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■ Abstract. We discuss, from a statistical point of view, some leading issues 
f-^ | that deal with the study of stellar populations in fully or partially unresolved 

, aggregates, like globular clusters and distant galaxies. A confident assessment 

■ of the effective number and luminosity of stellar contributors can provide, in 
this regard, a very useful interpretative tool to properly assess the observational 
bias coming from crowding conditions or surface brightness fluctuations. These 
arguments have led us to introduce a new concept of "photometric entropy" of 
a stellar population, whose impact on different astrophysical aspects of cluster 

. diagnostic has been reviewed here. 
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Galaxy surface brightness fluctuations and crowding effects in nearby star clus- 
ters are two related and well recognized features one has to deal with when 
Q H i observing partially or fully unresolved stellar systems. The latter effect turns 

out to be a severe problem, for instance, when probing the innermost regions of 
globular clusters in our own galaxy, and more generally in all those situations 
c/2 , when the stellar plot consists in fact of blended point sources, at different spatial 

scales. 

On the other hand, to a more detailed analysis, even the smooth sur- 
face brightness of distant galaxies actually is found to display some intrinsic 
"dumpiness", and this special property can be usefully exploited to derive in- 
for mation on their ful l y blen ded compos i ng st ellar populations, as first shown 
bv lTonrv &: Schneide r (1988) and lTonrvl (|l99ll b 

In this contribution, we would like to further extend the analysis of these 
two important issues, that deeply relates to overall characteristics of stellar 
aggregates and the way the latter are sampled by the observations. In particular, 
we will show that both effects are in fact a consequence of the same property 
of finiteness and discreteness of the composing stellar populations, and this will 
lead to a unified and more general definition of "photometric entropy" of a stellar 
population, a relevant concept whose impact on different astrophysical aspects 
of cluster diagnostic we want to review here. 

2. Observational constraints 

Our attempt to single out any "clean" information of a cluster stellar population 
crucially relies on how confidently we approach the ideal observing conditions, 
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Figure 1 . Upper panel - The ex pected trend of seeing variation with varying 
wavelength, after |Wynngj |l999|) . From empirical estimates and theoretical 
arguments, the FWHM of the seeing disk is found to scale approximately as 
FWHMscoing oc A -1 / 5 . Lower panel - Telescope resolving power with varying 
the aperture diameter, D, (right labels, in meters). As predicted by the 
classical Airy diffraction pattern, FWHMAiry oc X/D. Contrary to seeing, 
optics performance degrades when moving to long-wavelength photometric 
bands. 



that is those allowing us to fully resolve any star member in the system. Apart 
from the obvious constraint due to the physical distance of our target and its 
apparent size, there are at least three relevant external mechanisms that play a 
role to bias our results when carrying out imagery (or spectroscopy) of a stellar 
aggregate. 

• Photon noise directly relates to the signal-to-noise ratio (S/N) of our 
observations, the latter being constrained both by telescope size and, in 
case of ground-based observations, by sky brightness. In particular, we 
know that, for faint sources, S/N mainly scales with the telescope aperture 
(D), the exposure time (t exp ), and the sky surface brightness (ssb), in 
magnitudes per square arcsec, such as 

S/N oc (1/D) x (1/VW) x 10 a2ssb . (1) 

• Telescope resolving power is also an issue. It can be quantified by 
the full width at half maximum (FWHM) of the instrumental diffraction 
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pattern (the "Airy disk"), that is a measure of the minimum physical 
"spot" of a star on our imagery detector (or spectrograph slit), if we could 
observe outside the atmosphere. The Airy disk depends itself on D, but 
scales with the observing wavelength, as well, so that 

FWHM Airy ~ 4- ( 2 ) 



Finally, seeing conditions, in case of ground-based observations, drasti- 
cally depend on sky clarity and, opposite to the telescope resolving power, 
improve in general at longer wavelength, mainly due to a better spatial co- 
herence of the atmosphere convectio n layers with increasing A. Theoretical 
arguments and empirical estimates ( Wvnneiri999fi indicate that FWHM of 



the seeing disk scales as 

FWHM seeing a A" 1 / 5 . (3) 



While, to some extent, we could improve S/N of our observations (and 
"see more clearly") by arbitrarily increasing exposure time (although with some 
non- negligible technical limits due to detector saturation, electronic noise etc.), 
the other two effects are much more difficult (and expensive) to be optimized, 
as one would need either a bigger telescope and/or space-based observations to 
fully remove any atmosphere contamination. The combined action of telescope 
resolving power and seeing provides therefore the observational kernel that con- 
strains our ability to master the crowding conditions and fully resolve the stellar 
population (see Fig. P). 

For example, for the typical case of a ~ 50 pc globular cluster around the 
Andromeda galaxy (about 0.7 Mpc away) a mean angular separation of the 
order of 0.04" should be expected among the ~ 10 5 star members of the cluster. 
Such an equivalent resolving power would certainly not be achieved by means of 
any seeing-limited conventional telescope on the ground, rather demanding at 
least a 2.5 m orbiting telescope (as it is actually the case for the Hubble Space 
Telescope, indeed!). 



3. Effective numbers of luminous contributors and photometric en- 
tropy 

Although within the limits of the instrumental performance and seeing influence 
(both being, however, quantitatively assessed and, at least partially, recovered), 
one can still take advantage of the overall statistical properties of a stellar sample 
and gain valuable information, in any case, of its unresolved component. 

In this regard, photometric entropy adds an important tool to our analysis. 
As we will see in a moment, its definition has much to do with the concept of 
"effective number" of stellar contributors to the integrated luminosity of a stellar 
population. This qu antity has been first assessed in a series of previous papers 
( Buzzonilll989l. I1993T). providing the reference framework for further theoretical 
investigations (ICervifio. Luridiana h Castanderll20'ool : ICerviho et al l200li 12002: 
ICervifio fc Luridianall2004l ~ We will just sketch here some of the leading issues 
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Figure 2. The Poissonian distribution of the iVtot star cells in our "fair sta- 
tistical representation" of a stellar population. In average, each cell contains 
one star, according to the Poisson formula, p{x) — e _1 jx\. As a consequence, 
after repeated trials, the total number of stars, A to t, is expected to fluctuate 
by a{N tot ) = VAW 



of the theory and the main relationships among the relevant statistical quantities 
involved in our analysis. 

First of all, it is important to define operationally what we call a "fair 
statistical representation" of a stellar population. Ideally, we could assume to 
have a number N tot of cells, each to host, in average, one star to be supplied to 
the system through a stochastic process of Poissonian nature. After completion 
of N^t iterations, each cell will display a star distribution (rij) like in Fig. O 
and the whole system will contain, in average, J2f=i n j = ^tot stars. 

Due to the Poissonian distribution of nj, however, in repeated statisti- 
cal realizations of the population, the total number of stars will fluctuate by 

ff(iVtot) = 

If we assume, to a first analysis, that all stars have the same luminosity 
I, then the expected relative fluctuation of the global luminosity (Ltot) of the 
sample is: 



Ltot zZj nji IN 



l/VNot- (4) 



More likely, if lj is not a constant, and distributes according to a given 
luminosity function, then we could still retain eq. just replacing Ntot m the 
r.h. side of the formula with a more general quantity iV e g . The value of A^ e g can 
be regarded as an "effective" number of stellar contributors in the population; 
it could be demonstrated that, always, 



iVeff < Ntot- 



(5) 
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Figure 3. Logarithm of effective luminosity (£ c s) and sta r number (N P fr ) 
for a 15 Gyr SSP of solar mctallicity and Salpeter IMF, from lBuzzonil (|l993). 
Both quantities relates to the global luminosity (Ltot) of the stellar population 
through eq. H10|) . The model has been rescaled to a total SSP bolometric 
luminosity of 10 7 L©. 



Quite importantly, as the luminosity distribution of stars changes with wave- 
length, one has to expect that 

N cS = JVeff(A). (6) 

On a similar argument, the relative variance of the luminosity (Xtot) for the 
whole stellar population simply results: 



CT 2 (L tot ) _ Ej-n 



3<-j 



Ltot E, n 



4ff- (7) 



j '"3^3 



In the equation, -£ c ff is the "effective luminosity", that is basically a "mean" 
representative luminosity of the composing stars, using Wj = nj£j/L tot as a 
normalized weighting factor. 

It is interesting to note that eq. (J7|) is the key relation for the lTonrv Sz Schneider! 

( 1988f ) theory of galaxy surface brightness fluctuations. Actually, one major is- 



sue of the Tonry & Schneider method is that, relying on the observed galaxy 
flux /g a i, one can supply the empirical quantity 

^(/gal) = 4ff (8) 

/gal And 2 

(being d the galaxy distance), to be matched with the corresponding theoretical 
predictions for £ e g from population synthesis models, according to the distinctive 
properties of the (unresolved) galaxy stellar population. As a result, the method 
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Figure 4. A schematic representation of N e s, defined in our statistical ex- 
periment as the maximum number of "switched-on" cells, each containing one 
star of fixed luminosity £ e s, such as to provide the total observed luminosity 
of the population, Ltot, at a given photometric band. 

proves to be, in principle, a powerful distance indicator, leading to a direct 
measure of the galaxy distance modulus: 



(m — M) = 5 log d — 5 = 5 log 

By matching eqs. (j3J) and (JJJ), a relevant physical constraint should hold 
(see Fig. |3J) as: 

N cS x 4 ff = Ltot- (10) 

Accordingly, from a statistical point of view, iV e ff represents therefore the 
maximum number of bright stars, of constant luminosity t e &, allowed in a pop- 
ulation of -ZVtot members to provide the total luminosity L to t (see Fig. |1J). Ide- 
ally, this statistical definition assumes a "two-state" condition for sample stars, 
with "switched-on" objects of £ e s individual luminosity and the remaining 
(^Vtot — -Weff) "switched-off" stars with £ = 0. Alternatively, N c g can also be 
regarded as the maximum number of states (i.e. "switched-on" cells in our pre- 
vious example) available to the system, and this straightforwardly leads to a 
more general extension of the "entropy" (5) concept, that in our framework can 
now be defined as 

S = log(AW^tot)- (11) 
or, in its usual thermodynamical notation, 

AS = A log N cS (12) 

as we would better like to single out its variation within a stellar population 
rather than its absolute value. 

According to our normalization, as < N e s(X)/N to t < 1, S is always a 
negative quantity, depending on wavelength and tipping at zero, in case of a 
"maximum-entropy" stellar system, where all stars contribute with the same 
luminosity at a given photometric band. 
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Figure 5. The centermost 3° x 3° region of the An dromeda galaxy (M 31), 
as imaged by the ultraviolet space telescope GALEX l|Thilker et al.ll2005j) be- 
tween 1500 and 3000 A (left panel), and seen at optical wavelength (right 
panel). The image is courtesy of the GALEX scientific team, and made avail- 
able at the Project Web site http://www.galex.caltech.edu/ It is evident the 
increased "transparency" of the galaxy at ultraviolet wavelength, as a conse- 
quence of a drastically lower number of effective stellar contributors. See text 
for a discussion. 



4. Crowding effects and apparent opacity of stellar systems 

As shown in Fig. |SJ one striking feature of the iV e fj function (and its derived 
entropy quantity) is its wide change with wavelength. In par ticular, one see s 
that for a 15 Gyr simple stellar population (SSP) model, from iBuzzon I {1993), 
the entropy peaks around blue/visual wavebands (UBV), while the number of 
effective contributors dramatically drops by nearly five orders of magnitude when 
moving to short wavelength, below 3000 A. A similar (although much milder) 
trend is also evident at infrared wavelength, where N e Q smoothly decreases by 
about one dex, compared to the optical range. 

Scaled to a typical 10 11 L Q galaxy, this means that, by looking at mid- ul- 
traviolet wavelength, we have to expect a thin plot of some 10 4 UV-bright stars 
tracing the galaxy body, compared to a quite smooth surface brightness distribu- 
tion at visual wavelength, provided by about 10 10 effective stellar contributors. 
The Andromeda galaxy, as imaged in the v isible light and by the ultraviolet 
space telescope GALEX ((Thilker et al.l 12005), is a good example in this sense 
(see Fig. |5J) • 

The exact trend of iV e g (and S) vs. wavelength is a natural output of theo- 
retical codes for population synthesis, and it can easily be computed for a wide 
range of the distinctive evolutionary parameters for a stellar popula tion. De- 
tailed values for SSP grids of models can be found in iBuzzonj (Il993l. and W eb 
updates at http://www. bo. astro. it~eps/home. html ) and lCerviho et al.l (12002ft . 

The knowledge of the entropy level is of special importance in order to 
quantitatively assess the expected crowding conditions or rather the apparent 
"optical depth" , when observing a distant stellar system at different photometric 
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Figure 6. The apparent system opacity (r) in case of a SSP as in Fig. 
assuming to observe the stellar sample from ground and from space (up- 
per curves). Crowding conditions are mainly modulated by seeing, in case 
of ground-based observations with conventional telescopes, and by the tele- 
scope diffraction limit in case of space imaging. Note that calculations for 
ground-based observations have been stopped shortward of 3000 A given the 
atmosphere opacity. In addition, the two curves have been offset by an arbi- 
trary quantity in the logarithm domain to match at the V band. Due to an 
opposite trend of seeing and optics FWHM (see the two curves at the bottom, 
as labeled), the convolution with the apparent surface density of stars (N e s 
curve) leads to different performance outputs. In particular, ground-based 
observations sensibly improve when moving at longer wavelength, thanks to 
a more favorable seeing condition, while crowding conditions are expected to 
maintain roughly constant at the different photometric bands in case of space 
observations. The Johnson/Cousins photometric band system is reported at 
the top of the figure, as a reference. 



bands. The latter will depend in fact on the effective number of stellar contrib- 
utors convolved with the instrumental kernel (telescope diffraction pattern plus 
seeing PSF) that degrades our resolution leading to a smooth nearly "solid" 
surface brightness profile. 

In case of ground-based (i.e. seeing limited) imaging of a star cluster, for 
instance, the apparent optical depth (r) is expected to scale as 

iFWHM^ 

log r ground = S + log — — -f h const. (13) 

Cluster apparent area 

If, on the contrary, we are observing from space (that is fully exploting the 
telescope diffraction limit), then target opacity will scale as 

vrFWHMi irv 

log T space = S + log — — f h const. (14) 

Cluster apparent area 

In this framework, the difference between a "crowded field" and a "sur- 
face brightness fluctuation" is just a matter of apparent "opacity" of the stellar 
system, eventually depending whether r is respectively < 1 or > 1. 
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Figure 7. An illustrative example of surface brightness fluctuations in the 
globular clsuter M53 (alias NGC5024). The i? c -band frame has been taken 
at the 2.12m telescope of Cananea (Mexico) with a FOSC camera. Field 
is roughly 6.5' x 6.5' sampled with a 1" pixel size of on the sky. For our 
experiment, the cluster has been divided in four quadrants reporting the sky- 
subtracted integrated CCD counts for each partition, as labeled on the plot. 
See text for discussion. 

According to the trend envisaged in Fig. ^ the expected performance of 
our observations will therefore change with A as in Fig. El In general, we see 
that ground-based observations are more sensitive to the observing photomet- 
ric band as stars will more easily be picked up at longer wavelength, taking 
advantage of seeing improvement; on the contrary, our output will not be so 
critically constrained when observing from space, as the lower value of N e s at 
infrared wavelength barely compensates the poorer telescope resolving power 
(compared to the visible wavelength range), thus maintaining the crowding con- 
ditions roughly unchanged along the different photometric bands. 

5. Surface brightness fluctuations 

The previous arguments in our discussion make easier now a different statistical 
approach to the problem of surface brightness fluctuations in stellar systems. In 
particular, relying on our definition of "fair statistical realization" of a stellar 
sample (see Sec. 3), one could take advantag e of cluster and galaxy sym metry for 
a different and more general statement of the lTonrv &: Schneider! (|l988l ) method. 

For example, in case of a globular cluster, like the i? c -band CCD picture 
of M 53 shown in Fig. [7J one could assume that any quadrant of the system (by 
centering on the photometric barycenter of the cluster) collects a "statistically 
fair sample" of the whole cluster population, in force of the claimed central 
symmetry of the system. So, while, in average, each quadrant will display a 
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mean luminosity V = L to t/4, a scatter has to be expected for this quantity, 
according to l.h. side of eq. (J7|, simply measured as the variance of the four 
count determinations. In the specific case of Fig. [7|we have 



L' = 3138 counts 
a 2 (L') = 184 2 counts 



(15) 



The apparent (i? c -band) effective magnitude (m e ff) of the M53 stellar population 
becomes therefore 

m eS = -2.5 log = _ 2 .5 log + M R (M53), (16) 

being L tot = 4L' the apparent luminosity of the cluster (i.e. across the four 
quadrants, in CCD counts), to be calibrated in magnitude scale throu gh the 
known apparent i? c -band magnitude of M53 [M/j(M53) = 7.11, from iHarrisl 
(1996), dereddene d assuming a col o r exc ess E(B — V) = 0.02 and A(R C ) = 
2.6E(B — V) from lMoro &z Munaril ((2000t )]. Replacing the relevant quantities, 

for our cluster we have m e s = 14.78 mag; 

For the M53 stellar population, the lBuzzonil ( 19891 . H993T ) SSP models pre- 



dict and absolute (i? c -band) effective magnitude M e g- = —1.58 mag, assuming 
an age of ~ 15 Gyr f or the cluster, a metallicity [Fe/H] = —1.88 ± 0.2 dex 
( Santos fc Piattil 1200^ and a blue horizontal branch (HB) morphology. This 



eventually leads to an estimated distance modulus for M53 of 

m eS - M eff = 14.78 + 1.58 = +16.36 mag (17) 

in quite good a g reeme nt with the standard value of (to — M) = 16.25 mag 
from the lKirS) globular cluster catalog; this implies a ~ 5% relative 
uncertainty in the derived distance of the cluster with our method. 

5.1. Color fluctuations and stellar sampling 

Previous application of cluster "fair sampling" basically relies on the claimed 
symmetry of the system. In this regard, our adopted 4-quadrant partition is 
fully arbitrary, as any other geometrical combination, like slices of fixed angular 
aperture, or other fixed-size apertures symmetrically located around the center 
could work equally well in our method. 

As an interesting variant on this line, that may more suitably apply to fully 
unresolved galaxies, one could also consider to trace the luminosity variance 
along a given isophote of the surface brightness profile. This could easily be 
done, for instance, in case of ellipticals or face-on spirals. 

One important consequence of the finite luminosity sampled per pixel res- 
olution element, along a given galaxy isophote, is that some intrinsic variance 
of the apparent color has to be expected, as a consequence of a different pho- 
tometric entropy with varying the observing band. The amplitude of this color 
scatter depends on the amount of bolometric Lq sampled per pixel and the 
relative variation of iV e g at the different photometric bands. 

An example is done in Fig. |H1 where we computed the statistical variation 
of the integrated B-V color versus sampled luminosity. The B-V variance (in 
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Figure 8. The expected statistical scatter in the galaxy B-V color versus 
pixel luminosity sampling. The scatter amount ranges according to eq. (|19fl . 
High-S/N (i.e. S/N ^ 100) high-resolution galaxy imagery could in principle 
detect the effect. 



magnitude scale) derives from 

a 2 (B -V) = <j 2 (B) + <j 2 (V) - 2pa{B, V) (18) 

where a 2 (B) = cr 2 (V) = l/N^, and the covariance term -2pa(B,V) 

ranges between zero (if the B and V luminosity contributors can be assumed 
to be totally independent, that is for a correlation coefficient p = in eq. I18|) . 
and — 2a(B)a(y) if we assume a perfect positive correlation between the two 
photometric bands (i.e. assuming p = 1). As shown in Fig. |SJ the full variation 
range for the statistical fluctuation of the integrated color can eventually be 
written as 

| a(B) - a{V) \ <a(B-V)< \Ja 2 (B) + a 2 (V) (19) 

Quite interestingly, one sees from the figure that a detectable color scatter 
could be measured along a galaxy surface isophote providing to collect high- 
(S/N) (i.e. S/N S> 100) high-resolution imagery. For example, in case of a 
typical 10 11 L Q elliptical galaxy at the Virgo distance (d ~ 15 Mpc), one expects 
to sample roughly 10 5 L /px with a CCD of ~ 0.2" pixel size. This leads to a 
fully measurable a(B — V) scatter of some 0.01 mag (see Fig. |SJ). The possible 
cosmological relevance of this test is obvious, as from a measure of a color scatter 
in a distant galaxy we can derive an absolute estimate of the sampled luminosity, 
to be compared with the observed surface brightness and derive therefrom the 
luminosity distance. 



6. Diagnostic tools for high-resolution spectroscopy 

Following a substantially similar argument as for the color statistical scatter of 
previous section, one could also further expand the approach including spec- 
troscopy. Again, the change of 5 and N c q along wavelength is the key issue 
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Figure 9. Theoretical spectral energy distribution (upper plots) and effec- 
tive stellar contributors (N e g) (lower plots) for two SSPs of solar metallicity 
and Salpeter IMF according to Buzzoni (1989, 1993) population synthesis 
models. The models refer to a 2 Gyr SSP with red HB morphology, and a 
15 Gyr populatio n with blue HB, as labeled. T he new Uvblue theoretical 
library of spectra ijRodriguez-Merino et al.ll2005|) has been used to reproduce 
the spectral region around the striking Mgn and Mgi features, about 2800 A. 
In spite of a substantially indistiguishable difference between the spectra, the 
two SSPs display a large difference in terms of effective contributors and pho- 
tometric entropy, with a much higher statistical scatter in the spectral features 
expected for the older SSP. 



that leads us to expect some variance among the spectral features, even at close 
wavelength difference. 

A confident assessment of this effect could be extremely useful as an addi- 
tional interpretative tool to disentagle the well recognized degeneracy among the 
evolutionary parameters of stellar populations (see, e.g. the lo ng-lasting question 
of the "ag e- metallicity" dilemma as discussed, for instance, bv lRenzim & Buz zonil 
Il98ti and iBuzzonilllS . 

Figure is an illuminating example of this kind of problem. The two plots 
in the figure display the ultraviolet spectral energy distribution and effective 
number of stellar contributors for two SSPs of di fferent age and HB morphology. 
The models have been obtained by matching the lBuzzonil (1989) s ynthesis code 
with t he new Uvblue theoretical spectral library of iRodriguez- Merino et al.l 
( 2005T ). 
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In particular, a young (2 Gyr) population with red HB morphology is com- 
pared with an old (15 gyr) one, with a blue HB. One sees that, in spite of the 
striking age difference, the two spectra are nearly identical as, in both cases, 
we have a dominating component of warm stars (T e fj ~ 10 000 K) that emit in 
the ultraviolet range. However, while in the 2 Gyr SSP the UV-luminosity is 
provided by a large number of main sequence turn-off stars, in the 15 Gyr case 
we have a prevailing contribution from a few bright stars in the blue tail of the 
HB. 

This has a direct impact on the value of iV e fj, with a much larger scatter 
along wavelength for the oldest SSP. As a consequence, in the latter case we 
should likely expect a larger variance (about a factor of two) in the measure of 
any narrow-band spectrophotome tric index along the ga laxy isophotes [including 
bot h the popular Lick indices of IWorthev et all ( 1994 ) , in the optical range or 



the iFanelli et al.l (|l990l ) ultraviolet indices] , in a way very similar to what we 
have shown for the B-V color. 
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